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Abstract 

The paper contains integral representations for certain classes of 
exponentially growing solutions of second order periodic elliptic equa- 
tions. These representations are the analogs of those previously ob- 
tained by S. Agmon, S. Helgason, and other authors for solutions of the 
Helmholtz equation. When one restricts the class of solutions further, 
requiring their growth to be polynomial, one arrives to Liouville type 
theorems, which describe the structure and dimension of the spaces 
of such solutions. The Liouville type theorems previously proved by 
M. Avellaneda and F. -H. Lin, and J. Moser and M. Struwe for peri- 
odic second order elliptic equations in divergence form are significantly 
extended. Relations of these theorems with the analytic structure of 
the Fermi and Bloch surfaces are explained. 
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1 Introduction 



The topic of this paper stems from two sources. The first of them are 
representation theorems for certain classes of eigenfunctions of the Laplace 
operator in R", or equivalently, of solutions of the Helmholtz equation 

-An -k'^u = in M" , (1.1) 



where k € C* := C \ {0}. Such theorems for arbitrary solutions of ( 1.1 
were obtained in M? and in the hyperbolic plane by S. Helgason [22, p^ , 
and in by M. Hashizume et al. [^l|, M. Morimoto |37], and recently by 



S. Agmon We remark that it should also be possible to deduce similar 
results from the L. Ehrenpreis' fundamental principle. The zero set of the 
symbol of the operator in the left hand side of ( |1.1D is 

S = e C"|e^ = k^} , 

where = X]j=i?j- The L. Ehrenpreis' "fundamental principle" in the 
particular case of equation ( |1.1[ ) claims that any solution of ( |1.1D can be 
represented as a combination (i.e., an integral with respect to the parameter 
^) of the exponential solutions 

eg(x) := exp(i,^ ■ x), ^ G S, 

where ^ • x = X]j=i ^j^j (see the details and more precise formulation in [^] 
or [^]). The set S is an irreducible analytic subset of C", which is uniquely 
determined when A: ^ by its spherical subset 

Here S^~^ denotes the unit sphere in R"'. It is clear then that due to 
the uniqueness of analytic continuation, the exponential representation of 
solutions u{x) of (|1.1| ) should be reducible to one that involves the solutions 
e ^ with ^ G 5* only. Namely, consider the restriction mapping from functions 
analytic on the whole characteristic variety S to the sphere S. Due to the 
irreducibility and the uniqueness of analytic continuation, this mapping is 
one-to-one. Hence, there is a function space on the sphere S which is the 
isomorphic image of the space of all analytic functions on S. It follows that 
any hyperfunction (analytic functional) on S can be rewritten as a functional 
on S. Since the "fundamental principle" essentially expresses all solutions 



of (1.1) as applications of such analytic functionals to the analytic family of 



exponential solutions, we get our conclusion. 
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Now, depending on how fast the solution u{x) grows at infinity, the corre- 
sponding representing functional on S is actually a measure, a distribution, 
a hyper function, or a functional of a more general kind. For instance (see 
[Q, ^]), an arbitrary solution in M" can be represented as 

u{x)=<^{0,e^{x)>, (1.2) 

where (p{(,) is a functional on S which belongs to the dual space to the space 
£ := liuiR^oo £r{S^^^). Here for every i? > the Hilbert space £r{S^^^) 
is defined as follows: 

£r{S-'') := { V I V'(u;) = E"'.™f7]^^^™(^)' 

l,m ^ "2 ' 

\m£n--= (5^k™P)^/2<oo}, 
l,m 

where Yi^{uj) denote spherical harmonics, and (?(S""~^) is equipped with the 
inductive limit topology of lim/j^oo £r{S"'~^). 

The representation ( |1.2D can be formally rewritten as 

uix) = / e^{x)d(p{C). 
Js 

where is a suitable functional. 

The functional </) is a hyperfunction (analytic functional) on S if and 
only if for arbitrary e > the solution u{x) grows not faster than 

0{exp{{\lmk\ + e)\x\)) 

(see Q). One can also describe other classes of solutions, for instance, 
solutions which are represented by a distribution or a measure (see [§, 
36 1 and the references therein). 

As we have already mentioned, these results could be probably extracted 
from the "fundamental principle" | p!6| , pS| ]. The crucial factors are that S is 
sufficiently massive and S is irreducible, so S determines S uniquely. Be- 
sides, S is a rather simple analytic manifold. These features allow more or 
less explicit descriptions of the needed spaces of test functions and function- 
als. It is easy to understand that if S were reducible, it would not be possible 
to obtain the representation of all solutions using only ^ € S. The reason 
is that the solution with ^ that belongs to a component not touching S 
would not be representable this way. On the other hand, if one wants to 
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deal only with solutions growing not faster than 0{exp{\lmk\ + for 
all e > 0, then the irreducibility of S is not needed. In this case, it is only 
required that S is irreducible in a vicinity of S, so other components of S 
do not meet S. 

The "fundamental principle" was extended in ||3^ to solutions of certain 
growth (for instance, of exponential growth) of elliptic and hypoelliptic pe- 
riodic equations (see also the extensions of the results of |Q provided in 
[|40|). The role of the exponential solutions is played here by the so called 
Floquet-Bloch solutions (see Definition ||), and an analog of the character- 
istic manifold S is the variety F sometimes called the Fermi surface (see 
Definition ^). This raises the hope of finding representations similar to the 
ones discussed above for the more general case of a second order elliptic op- 
erator with periodic coefficients. This is not straightforward, however, due 
to several reasons. First of all, it is not that clear what should be a natural 
analog of S. An appropriate variety, as we explain later, is provided for sec- 
ond order equations by the analysis of the cone of positive solutions done by 



S. Agmon and by V. Lin and Y. Pinchover (see ||2|, 3C], and the references 
therein). The disadvantage is that the whole consideration must be done 
below the spectrum of the operator (more precisely, below the generalized 
principal eigenvalue Aq, see (|2.8| )). Secondly, proving the irreducibility of F 
happens to be a very hard nut to crack (this problem arises also in direct 



and inverse spectral problems, see for instance |y, |3l|, |3^). Fortu- 

nately enough, by appropriately restricting the growth of the solutions, one 
can sometimes work near a single irreducible component, and hence avoid 
proving the irreducibility of F. Consequently, we prove a representation the- 
orem (Theorem 18) that characterizes all the solutions which have integral 



expansion into positive Bloch solutions with a hyperfunction as a "measure" . 

The "fundamental principle" also suggests a point of view that is cru- 
cial for establishing representation theorems for solutions of equations with 
constant or periodic coefficients. Namely, it is to one's advantage to treat 
solutions of the original equation in the dual sense, i.e., as functionals on 
appropriate spaces that are orthogonal to the range of the dual operator. 
We adopt this approach throughout the paper. 

If one attempts now to further restrict the growth of solutions and con- 
siders the problem of the structure of all polynomially growing (or bounded) 
solutions, one arrives at the second topic of our study, Liouville type theo- 
rems. The classical Liouville theorem characterizes the space of all harmonic 
functions in of polynomial growth of order N. The validity of an ana- 
log of this classical theorem has been studied in many situations (see for 
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instance |14, 33, ^ for recent results, surveys, and further references). An 
interesting case was considered by M. Avellaneda and F. -H. Lin and 
also by J. Moser and M. Struwe |^8|. In these papers the authors dealt with 
polynomially growing solutions of a second order elliptic equation Lu = 
in divergence form with periodic coefficients and obtained a comprehen- 



sive answer (for related results see also [14, 33 1 and the references therein). 
More precisely, using the formalism of homogenization theory ||To| , |25|| , it 
was proved that any solution v of the equation Lu = in M" of polynomial 
growth is representable as a finite sum of the form 

= ^ x^pj{x), (1.3) 
j=(ii,- jn)eZ!^ 

where the functions Pj{x) are periodic with respect to the group of periods 
of the equation. Moreover, the space of all solutions of the equation Lu = 
of polynomial growth of order at most is of dimension hn,N, where 

fn + N\ f n + N-2\ 
hn,N:=[^ N )-[ N-2 ) (1-^) 

is the dimension of the space of all harmonic polynomials of degree < A^ in 
n variables. We will also use the notation 

(ln,N ■= [ N ) (1-5) 

for the dimension of the space of all polynomials of degree at most A^ in n 
variables. Notice that qn-i,N also coincides with the dimension of the space 
of all homogeneous polynomials of degree A^ in n variables, so in particular, 

hn,N = qn-l,N-l + Qn-l^N ■ 

We remark that the method of [^, |3^ can be slightly modified to provide 
an extension of this Liouville theorem for general second order elliptic equa- 
tions with periodic coefficients under the assumption that the generalized 
principal eigenvalue is zero (see Appendix |A| and also the recent paper |35], 
where a partial result of this type was independently obtained). 

One can make an observation that these Liouville theorems are actually 
of the same nature as the representation theorems discussed above. In this 
case the analog of the set S is the single point ^ = and the representing 
functional </> is a distribution supported at this point. 

Let us recall the following standard notion of Floquet theory (see [ p^ , 



30, 4 
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Definition 1 A solution u{x) representable as a finite sum of tlie form 



u(x)=e^'=-^ J2 x^Pj{x)\ (1.6) 

with functions Pj{x) periodic with respect to the group of periods of the 
equation is called a Floquet solution with a quasimomentum k E C". Here 

n 

k ■ X = hxi The maximum value of \j\ = ^ ji in the representation ( |1.6| ) 

1=1 

is said to be the order of the Floquet solution. Floquet solutions of zero 
order are called Block solutions. 



One sees that the representation ( |1.3D corresponds to a Floquet solution 
with a zero quasimomentum. A Liouville theorem of the type mentioned 
above implies in particular that the only real quasimomentum that can occur 
for the equation under consideration is A; = (modulo the action of the 
lattice reciprocal to the group of periods). We show in the present paper 
that the finiteness of the set of real quasimomenta for a periodic elliptic 
equation is equivalent to the finite dimensionality of the spaces of solutions 
having a given polynomial growth and to their representation similar to, 



albeit more general than (L3). This statement is very general and holds for 
any periodic elliptic equation (it is also true for hypoelliptic equations and 
systems, although we not address them here). If some additional information 
is available on the analytic behaviour of the dispersion relations, one can 
find the exact dimensions of these spaces (see Theorem E3) . We present in 



Theorem 28 some classes of second order equations (including Schrodinger, 
magnetic Schrodinger, and general second order elliptic equations with real 
periodic coefficients) for which one can achieve all these sharp results. We 
show that the problem of calculating the dimensions of the spaces of Floquet 
solutions of a given polynomial growth reduces to a purely function theoretic 
question and can be resolved in a very general setting (Theorem pX|). 

The proofs of the results of this paper are largely dependent upon the 
techniques of the Floquet theory developed in |3^] . 

The outline of the paper is as follows. The next section introduces nec- 
essary notations and preliminary results from the Floquet theory and the 
theory of positive solutions of periodic elliptic equations. In particular, we 
obtain a new general result (Theorem [lO| ) on the dimensions of the spaces of 
Floquet solutions, which plays crucial role in our approach to Liouville the- 
orems. Section |3| contains the proof of the integral representation (Theorem 
18) analogous to Theorem 5.1 in [Q]. In Section ^, we discuss Liouville type 
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theorems. In particular, Theorems and 28 are estabhshed. In order to 
make the reading of the paper easier, we postpone the proofs of all the tech- 
nical lemmas to Section ^. Some conclusions and remarks are provided in 
Section |^. The Appendix contains an alternative proof of a part of Theorem 
28 using the homogenization technique similar to the one used in ^ . 

Results of this paper related to Liouville theorems were presented in 
March 2000 at the University of Toronto and at the Weizmann Institute. 

When the paper was being prepared for submission, P. Li informed the 
authors that the statement of the third part of Theorem ^ for the special 
case of an operator of the form L = —'Yl, 0'ij{x)didj was simultaneously and 
independently obtained in |35| using homogenization formalism . 



2 Notations and preliminary results 

Consider a linear (scalar) elliptic partial differential operator P{x^ D) of 
order m in M*^, n > 2 (in some parts of the paper we will restrict the class 
of operators further). Here we employ the standard notation D = j-^. The 
ellipticity is understood in the sense of the nonvanishing of the principal 
symbol Pmix, °f the operator P for all ^ € \0. The dual operator (the 
formal adjoint) P* has similar properties. Here we use the duality provided 
by the bilinear (rather than the sesquilinear) form 

< 9,f >= J f{x)g{x)dx. 

We assume that the coefficients of P are smooth and periodic with re- 
spect to a lattice V in M". The smoothness condition can be significantly 
reduced (see the Section]^. In fact, so far we only need that both operators 
P and P* define Fredholm mappings between the Sobolev space and L2 
on the torus T" = M"/r. 

An additional condition is required that would guarantee the discreteness 
of the spectrum of the "shifted" operators -P(x, D + k) on the torus T" for 
all G C". We need to exclude the possible pathological situation when the 
spectrum of P on the torus coincides with the whole complex plane (like 
in the case of the operator exp(i (p) d/d4> on the circle). For instance, self- 
adjointness of P could be such a condition. Another example is second order 
uniformly elliptic operators of the form ( |2.5| ). For more sufficient conditions 
see for example |]l|]. 

In what follows, the particular choice of the lattice is irrelevant and can 
always be reduced to the case F = Z", which we will assume from now on. 
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We will always use the word "periodic" in the meaning of "F-periodic" . 

We denote by K = [0, 1]" the standard fundamental domain (the Wigner 
-Seitz cell) of the lattice T = Z", and hy B = [— tt, tt]"' the first Brillouin 
zone, which is a fundamental domain of the reciprocal (dual) lattice T* = 
(27rZ)". We identify F-periodic functions, in the natural way, with functions 
on T". 

We introduce now the set that plays the role of the characteristic variety 
E discussed in the introduction. 

Definition 2 The (complex) Fermi surface Fp of the operator P (at the 
zero energy level) consists of all vectors G C" (called quasimomentd) such 
that the equation Pu = has a nonzero Bloch solution u{x) = e^^'^p{x), 
where p{x) is a F-periodic function. 

It would be useful later on to realize that in this definition the positivity 
of the solution is not required, and in fact the solution is usually complex. 

In many cases, it is convenient to introduce a spectral parameter A. This 
leads to the notion of the Bloch variety: 

Definition 3 The (complex) Bloch variety Bp of the operator P consists 
of all pairs (A;, A) S C"^^ such that the equation Pu = Xu has a nonzero 
Bloch solution u(x) = e^^'^p{x) with the quasimomentum k. 

It is clear that the Fermi surface is just the projection onto the /c-space 
of the intersection of the Bloch variety with the hyperspace A = 0. 

One can consider the Bloch variety Bp as the graph of a (multivalued) 
function X{k), which is usually called the dispersion relation. Then the 
Fermi surfaces become the level surfaces of the dispersion relation. Since 
the spectra of all operators P{x, D + k) on the torus T" are discrete, we can 
single out continuous branches Aj of this multivalued dispersion relation. 
These branches are usually called the band functions (see Q ) . 

The following analyticity property of the Fermi and Bloch varieties is 
important: 

Lemma 4 p^ , Theorems 3.1.7 and The Fermi and Bloch varieties 

are the sets of all zeros of entire functions of a finite order inC^ and C"+^, 
respectively. 

Another property of the Bloch and Floquet varieties that we will need 
later is the relation between the corresponding varieties of the operators P 
and P*. 
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Lemma 5 Theorem 3.1.5] A quasimomentum k belongs to Fp* if and 
only if —k G Fp. Analogously, {k,X) € Bp* if and only if (— fc, A) € Bp. 
In other words, the dispersion relations X{k) and X*{k) for the operators P 
and P* are related as follows: 

\*{k) = X{-k). (2.1) 

We note that the Fermi surface Fp is periodic with respect to the re- 
ciprocal lattice r* = (27rZ)". Therefore, it is sometimes useful to fac- 
tor out the periodicity by considering the (analytic) exponential mapping 
p : C" ^ (C*)", where 

z = p{k) = p{ki, ... ,kn) = (expifci, . . . ,expzfc„). 

This mapping can be identified in a natural sense with the quotient map 
^ C"/r*. We also introduce the complex torus 

r = p(R'^) = {zeC"I \zj\ = l,j = l,2,...,n}. (2.2) 

Definition 6 The image <I>p = p{Fp) of the Fermi surface Fp under the 
mapping p is called the Floquet surface of the operator P. 

The reader familiar with the Floquet theory immediately recognizes the 
Floquet surface as the set of all Floquet multipliers of the equation Pu = 0. 

The main tool in the Floquet theory is an analog of the Fourier transform 
(see js^, Section 2.2], ||4^), which we will call the Floquet transform Li (it is 
sometimes also called the Gelfand transform): 

f{x)^Uf{z,x) = Y,f{^-l)^\ zeiCT, (2.3) 

where we denote z'^ = zj^ z^ ■ ■ ■ zj^ . 

It is often convenient to use for the Floquet transform U the quasimo- 
mentum coordinate k instead of the multiplier z, where 

z = p{k) = (exp iki , . . . , exp ik^) • 

We need to recall now some definitions from [^]. For a point z G 
(C*)", we denote by Em,z the closed subspace of the Sobolev space H^{K) 
formed by the restrictions of functions v G H^J^^) that satisfy the Floquet 
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condition v{x + 7) = z'^v{x) for any 7 G F. One can show (see Theorem 
2.2.1 in @) that 



forms a holomorphic sub-bundle of the trivial bundle (C*)" x H"^{K). As 
any infinite dimensional analytic Hilbert bundle over a Stein domain, it 
is trivializable (see Theorems 1.3.2, 1.3.3, and 1.5.23 in |3^). One can also 
notice that for m = the bundle Sq coincides with the whole (C *)" x L'^{K). 
The following standard auxiliary result for the transform lA collects sev- 



eral statements from Theorem XIII. 97 in [^] and Theorem 2.2.2 in |3C]: 

Lemma 7 1. For any nonnegative integer m the operator 

U : F"(M") ^ L\T,£m) 

is an isometric isomorphism, where Lp'{T,£rn) denotes the space of 
square integrahle sections over the complex torus T of the bundle £m, 
equipped with the natural topology of a Hilbert space. 

2. Let the space 

9"^ = {/ G iiZSR'')\sup\\f\\nrr^^K^^^ exp(6|7|) < ^ , V6 > 0} 
7Gr 

be equipped with the natural Frechet topology. Then 

z^:e™^r((c*r,f„) 

is an isomorphism, where r((C*)",£'m) is the space of all analytic 
sections over (C*)" of the bundle Sm, equipped with the topology of 
uniform convergence on compacta. 

3. Let the elliptic operator P be of order m. Then under the transform 
lA the operator 

P : H'"\mP) ^ L^iW") 

becomes the operator of multiplication by a holomorphic Fredholm mor- 
phism P{z) between the fiber bundles Sm o,nd 8q. Here P{z) acts on 
the fiber of Sm over the point z G T as the restriction to this fiber of 
the operator P acting between H"^(K) and L'^{K). 
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Here is another standard way of looking at the morphism P{z). Let 
z = exp z/c, then commuting with the exponent exp ik ■ x one can (locally) 
trivialize the bundle 8m reducing it to the trivial bundle with the fiber 
H^{T'^), where as before T" = R"/r. At the same time the operator P{z) 
takes the form P{x,D + k) between Sobolev spaces on the torus T". 

Let us discuss the structure of the Floquet solutions (see Definition |l|) 
and of functions of Floquet type (|1.6| ) in general. For illustration, consider 
the constant coefficient case, where the role of the Floquet solutions is played 
by the exponential polynomials 



\j\<N 



It is well known that, considered as distributions, all such functions are 
Fourier transformed into distributions supported at the point {—k). More- 
over, the converse statement is also true. A simple but extremely important 
and relatively unnoticed observation is that under the Floquet transform, 
each Floquet type function of the form ( [L.6| ) corresponds, in a similar way, 
to a (vector valued) distribution supported at the quasimomentum (— fc). 
We collect below this fact and some other previously known properties of 
Floquet solutions, as well as a new result on the dimensions of the spaces of 
such solutions, which will play the crucial role in establishing the Liouville 
type theorems. 

First of all, every Floquet type function n (see (1.6)), being of exponential 



growth, determines a (continuous linear) functional on the space If, 
additionally, it satisfies the equation Pu = for a periodic elliptic operator 
of order m, then as such a functional it is clearly orthogonal to the range of 
the dual operator P* : 0™ — > 0*^. According to Lemma 0, after the Floquet 
transform any such functional becomes a functional on F ((C *)" , fo), which 
is orthogonal to the range of the Fredholm morphism P*{z) : 8m 8q 
generated by the dual operator P* : 0™ Q^. We are now ready to 
formulate the following auxiliary result. 



Lemma 8 1. A continuous linear functional u on 0" is generated by 



a 



function of the Floquet form (1^) with a quasimomentum k if and 



only if after the Floquet transform it corresponds to a functional on 
r((C*)",£^o) which is a distribution cp that is supported at the point 
V = exp(— ifc), i.e. has the form 



lil<JV \ 



/er((c*)",fo), 
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where qj € L?'{K). The orders N of the Floquet function (1-t) and of 
the corresponding distribution (p are the same. 



2. Let ttk be the dimension of the kernel of the operator 

P{x, D + k): //'"(T") ^ L2(T"). 

Then the dimension of the space of Floquet solutions of the equation 
Pu = of order at most N with a quasimomentum k is finite and does 
not exceed akqn,N- 

The estimate on the dimension given in the second part of Lemma ^ is 
very crude and in many cases can be significantly improved. In fact, as the 
following theorem shows, we obtain an explicit formula for the dimension 
of the space of Floquet solutions with a given quasimomentum in the case 
of a simple eigenvalue. This theorem seems to be new and constitutes the 
crucial part of the Liouville theorem proved in Section |^ (Theorem 23). 



In order to formulate this result, we need to prepare some notions and 
notations. 

Definition 9 Let Q be a homogeneous polynomial in n complex variables. 
A polynomial p{x) in M" is called Q -harmonic, if it satisfies the differential 
equation Q{D)p = 0. 

Let V denote the vector space of all polynomials in n variables, and 
let Pi be the subspace of all homogeneous polynomials of degree /. Denote 

N 

by Vn = ® Pi the subspace of all polynomials of degree at most N . So, 

1=0 

oo 

V = @ Pi- If Q{k) is a nonzero homogeneous polynomial of degree s, then 

1=0 

the differential operator Q{D) : Pi^g Pi is surjective for any / (this simple 
statement will also follow from the proof of the theorem below). Hence, the 
mapping Q{D) : V ^ V has a (nonuniquely defined) linear right inverse R 
that preserves the homogeneity of polynomials. 

Theorem 10 Assume that zero is an eigenvalue of algebraic multiplicity 1 
of the operator P{x, D + ko): //"'(T") L'^{T') on the torus T". Let \{k) 
be an analytic function in a neighborhood of kQ such that X{k) is a simple 
eigenvalue of the operator P{x, D + k): i7'"(T") ^ ^^(T") and A(feo) = 0. 
Let 

oo 

\{k) =Y,Hk - ko) 
l=lo 
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be the Taylor expansion of X{k) around the point ko into homogeneous poly- 
nomials such that is the first nonzero term of this expansion. Then for 
any G N the dimension of the space of Floquet solutions of the equa- 
tion Pu = in of order at most N and with the quasimomentum ko is 
equal to the dimension of the space of all Xi^ -harmonic polynomials of de- 
gree of at most N. Moreover, given a linear right inverse R of the mapping 
Xig{D) : V ^ V that preserves homogeneity, one can construct an explicit 
isomorphism between these spaces. 

Proof. It is sufficient to consider the case k^ = 0, since the general case 
reduces to this by a change of variables. Consider the operator family 

A{k) = P*{x,D- k) - X{-k) : i7™(T") ^ L^{T'') 

which is analytic in a neighborhood of 0. At each point k of this neighbor- 
hood A(k) has by the construction a one-dimensional kernel. Then, accord- 
ing to Theorem 1.6.13 in there exists an analytic non-vanishing vector 
V'(fe) G KerA{k). In other words, P*{x,D — k)ijj{x,k) = X{—k)ip{x, k). Let 
us choose a closed complementary subspace M to KerA{0) in H"^{T^). 
Then it is complementary to Ker A{k) in a neighborhood of 0. Since 
P*{x, D) has zero kernel on M and is Fredholm, we conclude that P*{x,D — 
k) has zero kernel on M for all A; in a neighborhood of 0. We denote by Il{k) 
the closed subspace in L^(T") defined as Il{k) = P*{x,D — k){M). Applying 
Theorem 1.6.13 of |3^ again, we conclude that n(A:) depends holomorphi- 
cally on A: in a neighborhood of (i.e., forms a Banach bundle) and hence it 
is complementary to Ker A{k). Representing now the operator P*{x,D — k) 
in the block form according to the decompositions 



and 



we get 

P*{x,D-k) 



fjm^jn^ = M ® Ker A{k) 

L2(T") =Il{k)®KerA{k), 
B{k) 



X{-k) 



where B{k) is an analytic invertible operator-function between M and I\.{k). 
If we now have a functional (j) on r((C*)"',fo) supported at t; = exp(O), 
such that it is orthogonal to the range of the operator of multiplication by 
P*{k), then it must be equal to zero on all sections of the bundle n(A;). 
This means that the restriction of such functionals to the sections of the 
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one-dimensional bundle KerA{k) is an one-to-one mapping. This reduces 
the problem to the following scalar one: find the dimension of the space 
of all distributions of order N supported at the origin such that they are 
orthogonal to the ideal generated by A(— fc) in the ring of germs of analytic 
functions. One can change variables to eliminate the minus sign in front of 
k. Due to the finiteness of the order of the distribution, the problem further 
reduces to the following: find the dimension of the cokernel of the mapping 

Ajv : :Pjv ^ :Piv , 

where Ajv(p) for p € is the Taylor polynomial of order at of the 
product X{k)p{k). Let us write the block matrix Aij of the operator A^v that 

N 

corresponds to the decomposition Vn = ® Pi- It is obvious that Aij = 

1=0 

for i — j <Iq. For i — j >Iq the entry Ajj is the operator of multiplication 
by acting from Pj into Pi. Since A^q ^ 0, for i — j = Iq the operator 
Aij of multiplication by Xig has zero kernel. Being interested in the cokernel 
of Ajv, we need to find the kernel of the adjoint matrix A^. The adjoint 

N 

matrix acts in the space P^*, where P* can be naturally identified with 

1=0 

the space of linear combinations of the derivatives of order I of the Dirac's 
delta-function at the origin. Here we have A*j = for j — i < Iq, and for 
j — i > lo the entry A^j is the dual to the operator of multiplication by 
Xj-i acting from Pj into Pj. In particular, since for j — i = Iq the latter 
operator is injective, we conclude that the operators A*^ arc surjcctivc. This 
enables one to find the dimension of the kernel of the matrix A|^ and even 
to describe its structure. Namely, let 

TV 

V' = (^o,...,V'iv)G0pr 

1=0 

be such that A*'tp = 0. Due to the triangular structure of A^, it is easy to 
solve this system. Indeed, it can be written as follows: 

J2 A*j^l;j = 0,i = 0,...,N-lo. 
j>i+lo 

Taking the Fourier transform, we can rewrite this system in the form 

Xj^i{D)^j = 0,i = 0,...,N -lo, 

j>i+lo 
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where ip denotes the Fourier transform of ip. Therefore, ^lJj is a homogeneous 
polynomial of degree j in M". For i = N — Iq we have 

This equality means that ipN can be chosen as an arbitrary A^g-harmonic 
homogeneous polynomial of order A^. Moving to the previous equation, we 
analogously obtain 

Az3(D)^^ + Xi,+iiD)^ = 0, 
or _ 

The right hand side is already determined, and the nonhomogeneous equa- 
tion, as we concluded before, always has a solution, for instance 

-i?(Ai„+ip)^) . 

This means that 

is a A/g -harmonic homogeneous polynomial of order — 1. We see that the 
solution ^pN-l exists and is determined up to an addition of any homoge- 
neous A;y-harmonic polynomial of degree — 1. Continuing this process 
until we reach i/jq, we conclude that the mapping 

where 

i>j 

establishes an isomorphism between the cokernel of the mapping A^r and 
the space of A/g-harmonic polynomials of degree at most A^. This proves the 
theorem. □ 

In the cases of the simplest structures of the Taylor series, the theorem 
implies the following: 

Corollary 11 Under the hypotheses of Theorem the following hold: 

1. If kQ is a noncritical point of the band function X{k), then the dimen- 
sion of the space of Floquet solutions of the equation Pu = in M" of 
order at most N with a quasimomentum is equal to the dimension 
Qn-i,N of the space of all polynomials of degree at most N in M"""*^. 
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2. If the Taylor expansion of the band function X{k) at a point starts 
with a nondegenerate quadratic form, then the dimension of the space 
of Floquet solutions of the equation Pu = in M"' of order at most N 
with a quasimomentum is equal to the dimension hn,N of the space 
of harmonic (in the standard sense) polynomials of degree at most N in 
M". In particular, this condition is satisfied at nondegenerate extrema. 

In both cases an isomorphism can be provided explicitly as in the previous 
theorem. 

Proof: 1. By our assumptions, the Taylor expansion of A(A;) starts with a 
nonzero hnear term Ai(A;) = X]j=i^j%) ^ The corresponding differ- 
ential operator is 

Q d ^ d 

where aj and /3j are real. Consider first the case when the vectors a = 
(oLj) and /? = are collinear. Then Ai(D) becomes ToZ^Tj^, where 
7o ^ is a complex number and 7 = (7^) is a nonzero real vector. A linear 
change of coordinate system brings Ai(D) to the operator (up to an 
irrelevant constant factor). Thus, the Ai-harmonic polynomials are exactly 
those independent on x\. Invoking Theorem 0, we get our conclusion in this 
case. Consider now the situation when a and /3 are linearly independent. 
Then a linear change of variables brings Ai to the form d/cfz, where z = 
Xi + ix2- Since any polynomial in variables (xi, ...,Xn) is a polynomial of the 
same degree in (z, z, X3, j;„), the Ai-harmonic polynomials are the ones 
depending on (z, X3, a;„) only (i.e. the ones analytic in z). This again 
reduces the number of variables to n — 1. 

2. By our assumptions, the first nonzero homogeneous term is a nondegen- 
erate quadratic form X2{k — ko), which is reducible to the sum of squares 
of coordinates by a linear change of variables. Therefore, in the new co- 



ordinates X2{D) = —A. Using Theorem 10, we obtain the desired result. 

□ 

In the remaining part of this section we restrict further the form of the 
operator. Namely, we consider now second order operators with real periodic 
coefficients of the form 

n n 

L = - ^ aij{x)didj + ^ bi{x)di + c(x). (2.5) 

i,j=l 4=1 
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It is assumed that the uniform elhpticity condition 



n 

2 



^ aij{x)CiCj >a^C 

i,j=l i=l 

is satisfied for all x, C £ IK", where a is a positive constant. 

For such operators, we introduce the function that will play the crucial 
role in our considerations. Its properties were studied in detail in [0], fSf 



and |41|. Consider the function A(^) : M" M defined by the condition 
that the equation 

Lu = A{(,)u 
has a positive Bloch solution of the form 

u^{x) = e^-^p^(x), (2.6) 

where p^{x) is F-periodic. 

Lemma 12 1. The value A(^) is uniquely determined for any ^ G M". 

2. The function A(^) is bounded from above, strictly concave, analytic, 
and has a nonzero gradient at all points except at its maximum point. 

3. Consider the operator 

L{^) = e-^ ^'Le^ '' = L{x, D - iO 

on the torus T" . Then A(^) is the principal eigenvalue of L[£^) with a 
positive eig enf unction p^. Moreover, A(^) is algebraically simple. 

4. The Hessian o/A(^) is nondegenerate at all points. 

One should note that since the function A(^) is analytic, it is actually 
defined in a neighborhood of M" in C". This remark will be used in what 
follows. 

Let us denote 

Ao = maxA(n. (2.7) 

It follows from ||2|, |3^ that an alternative definition of Aq is 

Ao = sup{A G M I a-u > such that (L - A)m = in R"}, (2.8) 



17 



and that in the self-adjoint case Aq coincides with the bottom of the spec- 
trum of the operator L. The common name for Aq is the generalized principal 
eigenvalue of the operator L in M". 

We win often need to assume that Aq is either nonnegative or strictly 
positive. In the self-adjoint case such an assumption has a clear spectral 
interpretation. In the next lemma, we provide some known sufficient con- 
ditions for the nonnegativity or positivity of Aq for operators of the form 



Lemma 13 Consider an operator L of the form 

1. Aq > if and only if the operator L admits a positive (super) solution. 
This condition is satisfied in particular when c{x) > 0. 

2. Aq > if and only if the operator L admits a positive solution of the 
form ^Ij). 

3. Aq = if and only if the equation Lu = admits exactly one normal- 
ized positive solution in M". 

4- If c{x) = 0, then Aq = if and only if J b{x)ip{x) dx = 0, where tp 

is the principal eigenfunction of L* on T" . In particular, divergence 
form operators satisfy this condition. 

5. Let ^ G M", and assume that u^{x) = e^'^p^{x) and u*i^ are positive 
Bloch solutions of the equations Lu = and L*u = 0, respectively. 
Denote by ip the periodic function u^u*_^ . Consider the function 

n 

bi{x) = bi{x) -2^aij{x){^j + {p^{x)y^djPi:{x)}, 
i=i 

and denote 

7 = (7i,-- - >7n) := ( j bi{x)ij{x)dx,. . . , j 6„(x)V'(x) dx). 

Then Aq = if and only if j = 0. 

Let us discuss also some additional properties that will play an important 
role in the sequel. Assume that Aq > 0. Then Lemma 12 implies that the 
zero level set 

H = {e G M"| A(e) = 0} (2.9) 
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is either a strictly convex compact analytic surface in of dimension n — 1 
(this is the case if and only if Aq > 0) , or a singleton (this is the case if and 
only if Aq = 0) . The manifold H consists of all ^ G such that the equation 
Lu = admits a positive Bloch solution u^{x) = e^'^p^{x). Moreover, the 
set of all such positive Bloch solutions is the set of all minimal positive 
solutions of the equation Lu = in M" [|2|, 3£]. It is also established that a 



function u is a positive solution of the equation Lu = in M" if and only if 
there exists a positive finite measure /u on H such that 

uix) = lM^)MO- 

o 

We denote by G the convex hull of H, and by G its interior. Note that 

o 

if Aq > then Aq = if and only if H = G and hence G= 0- 

Lemma 14 Suppose that Aq > 0. There exists a neighborhood W of G in 
C" and an analytic function 

W3(^p^{-) eH^iJ'') 

such that for any ^ E the function of x 

U(^{x) = exp(^ • x)p^{x) 

is a nonzero Bloch solution of the equation Lu = A(^)n with a quasimomen- 
tum —i£,. Moreover, one can choose the function p in such a way that it is 
positive for all ^ G H. 

Comparing the definitions of H and of the Fermi surface F^, it follows 
that 

-zH C Fl. 

The next lemma specifies further the relation between these two varieties: 

Lemma 15 Let Aq > 0. Then 

1. The intersection of the complex Fermi surface Fl with the tube 

T = {k eC\ Imk = {Imki,...,Imkn) G -G} (2.10) 

coincides with the union of the surface — iH with its translations by 
the vectors of the reciprocal lattice T* , i.e. consists of vectors k = 
—i^, + 7 where ^ G H and 7 G T*. Moreover, up to a multiplicative 
constant, any nonzero Bloch solution with a quasimomentum in the 
above intersection is a positive Bloch solution. 
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2. If Aq > 0, then the intersection of Fl with a sufficiently small neigh- 
borhood of —ir. is a (smooth) analytic manifold that coincides with the 
set of zeros of the function A{ik). 

Analogously to the definition of the Floquet surface ^ = ^l, we define 
the surface 

^ = p{-iE) = {z\z = {exp^i,... ,expCn), ^ € H} , (2.11) 
and the tubular domain 

V = piT), (2.12) 

where T was defined in (|2.10 ). The results of Lemmas and |l5| can be 
rephrased in terms of these objects: 

Lemma 16 Let Aq > 0. Then 

1. $ny = 

// Ao > 0, then 

2. The intersection of ^ with a sufficiently small neighborhood of ^ is a 
(smooth) connected analytic manifold. 

3. The intersections of ^ with neighborhoods of the tube V form a basis 
of neighborhoods of ^ in ^. 

4- For a sufficiently small neighborhood of ^ in ^ there exists an 
analytic function p : $e — > //^(T") such that for any z £ the 
function of x 

Uz{x) = z^p{z, x) 
is a nonzero Bloch solution of the equation Lu = 0. 
We will also employ the following lemma: 
Lemma 17 Consider an operator L of the form I{2.d0 

1. Assume that c{x) > 0. Then the only solutions of the equation Lu = 
of the type exp(ik ■ x)p{x), where k G and p is a T -periodic function 
are the constants. If such a nontrivial solution exists, then c{x) = 0, 
and A(0) =0 (i.e. G 

2. Suppose that the operator L admits a positive periodic supersolution 
ip € C^'"(M"). Assume that v{x) = exp(iA; • x)p{x) is a nontrivial 
solution of the equation Lu = 0, where k G M" and p is a periodic 
function. Then there exists C G C such that v = Cip, the function ip 
is a positive periodic solution, and A(0) =0. 
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3 Representation of solutions by hyperfunctions 



The main result of this section (Theorem |l8|) is analogous to Theorem 5.1 
in which characterizes the class of solutions of the Helmholtz equation 
that can be represented by means of hyperfunctions on S (see also the in- 
troduction of our paper). 

In order to state it, we need to introduce a new object. Let us denote 
by h{Lj), u) € S'^~^ the indicator function of the convex set G. Namely, 

h{Lo) =sup {uj ■ (3-1) 

where w • ^ = Yl^=i ^j^j inner product in M". The next Theorem will 

be stated in terms of this function. 



Theorem 18 Suppose that Aq > 0. Let u be a solution of the equation 
Lu = in satisfying for any e > the estimate 

\u{x)\ < Ce exp ((/i(x/ |x|) + e) \x\) , (3.2) 

where C ^ is a constant depending only on e and u. Then u can be represented 
as 



u{x) =< /i(^),U5(x) >, 



(3.3) 



where uc is the analytic positive Bloch solution corresponding to ^ G 



(see 
The 



Lemma I4), and is a hyperfunction (analytic functional) on [ 
converse statement is also true: for any hyperfunction /x on E, the function 
u{x) in (\3.3^ is a solution of the equation Lu = in M" which satisfies the 



growth condition {3.i). 



Remark 19 Using a standard elliptic argument it follows that the pointwise 
growth condition (^]^) is equivalent to the growth condition 



u{x) exp (— {h{x/ \x\) + e) \x\) e L^ 



(3.4) 



Proof of Theorem 1^ : Assume first that a solution u has the representa- 
tion ( |3.3D . We need to prove that u satisfies the growth condition (3.2). Due 
to the real analyticity of u ^ with respect to ^ and according to lemmas ij and 
15, can be extended to an analytic vector function u^{x) = exp{^-x)p^{x) 
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on an e- neighborhood Ue of H in iFj^. Since is a hyperfunction (analytic 
functional) on H, we have an estimate 



Hence we have 



\u{x)\ < C £ ™ax |iig(x)| 



\uix)\ < Ceuiax e^-^ = C',el^l('^(^/l^l)+^(^)) , 



where lime__*o <^(^) = 0, which gives (|3.2D . 

Suppose now that u satisfies ( |3.2D . We need to prove that u can be 
represented as in ( |3.3[ ). Let Ge be the e- neighborhood of G and he = h + e 
be the indicator function of G^. Consider the following Frechet spaces of 
test functions: 

Wm,e = HJ^^MI I < >m,5 < OO, V < <5 < e} , 

where 

< (t> := sup { I mH^^K+j) e^'^^^^/I^DI^I) } . 

It is obvious that the operator L* maps continuously W2,e into Wo,e- Con- 
sequently, the linear functional 

<u,(f>>:= / u{x)4>{x)dx 

is continuous on the space Wq^s for any s > 0. Since Lu = 0, Schauder 
elliptic estimates together with the periodicity of the operator show that 
estimates similar to (|3.2| ) hold also for the derivatives of u. One observes 
that n is a continuous functional on Wo,e which annihilates the range of the 
operator L* : W2,e — > W^o.e- Now Floquet theory arguments analogous to 
the ones used in Section 3.2] can be applied to yield ( |3.3| ). Let us make 
this part more precise. 

Our first goal is to obtain a Paley- Wiener type theorem for the Floquet 
transform in the spaces Wm,e- Let us denote by the domain in (C*)" 

Ve = {z = (zi, . . . , Zn) G (C *)" I Zj = exp ikj such that 

Im k = {Im ki, . . . , Im kn) G ( — Gg) } . 

and let 

V: = {Z={Z,,... ,Zn)G {C*r I = [Z^\ ...,Z-')e Ve} , 

The domains V,. form a basis of neighborhoods of the tube V, where V is 
defined by ( p. 12 ). The following statement is a Paley- Wiener type theorem 
for the transform U which is suitable for our purpose. 
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Lemma 20 1. The operator 

U:Wm,e^r{V:,£m) 

is an isomorphism, where T{V*,£m) is the space of holomorphic sec- 
tions over V* of the bundle £m, equipped with the topology of uniform 
convergence on compacta. 

2. Under the transform U, the operator 

L* : W2,e ^ Wo,e 

becomes the operator C{z) of multiplication by a holomorphic Fredholm 
morphism between the fiber bundles £2 and £q: 

nvt,£2)''^^ nvt,£o)- 

Here C{z) acts on each fiber of £2 as the restriction to this fiber of the 
operator L* acting between H'^{K) and L^{K). 

Let us choose a value eo > such that the intersection of $ with Ve is 
smooth and connected. This is possible according to Lemma |l^. From now 
on, we will only consider the values < e < Eq. 

Since the image Uu of the solution u under the Floquet transform U is 
a continuous linear functional on r(V"*,£'o) which is in the cokernel of the 
operator 

r(y;,^2) -^r(F*,^o), 

our task is to describe all such functionals. Several theorems of this kind 



were proven in [30|. In our current situation such a representation can be 
obtained rather easily, due to the simplicity of the structure of the Floquet 
variety inside V^- Namely, let Uz{-) = z^p{z, •) be the Bloch solution of the 
equation Lu = introduced in Lemma |l^. Let also 7^($£) be the space 
of holomorphic functions on = ^ DVe equipped with the topology of 
uniform convergence on compacta. We introduce the mapping 

t:T{V:,£o)^n{^e) 
which for any section f{z,x) of the bundle £0 produces 

tf{z) =< f{z^^,-),Uz >= / f{z~^,x)uz{x)dx. 



Here z ^ = (z/, . . . , z„ ^] 
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Lemma 21 Let < e < Eq, where Eq is the value defined above. Then 
the mapping t is a topological homomorphism and the following sequence is 
exact: 

T{V:,£2) T{V*,£o) ^ n{^e) - 0. 

This lemma practically finishes the proof of the theorem. Namely, the 
solution u after the Floquet transform leads to a continuous linear functional 
on r(V^*,£'o) that annihilates the range of the operator of multiplication 
by C{z). Lemma ^ implies that such a functional can be pushed down 
to the space 7Y(<I>£). Since this functional, due to the estimate (|3.2|), is 
continuously extendable to 7i{^ for arbitrarily small values of e, it is in 
fact a hyperfunction (analytic functional) /U on $ = f] <!>£. Hence, the 

e>0 

action < n, > of the functional n on a function (j) G Wo.e can be obtained 
as 

<u,(t) >=< fj,{z),t{z){U4>) > . 

Applying now the explicit formulas for the transforms U and t, one arrives 
to the representation (|3.3|). Indeed, 



t(U<t>){z) = j U(t){z ^,x)uz{x)dx (3.5) 

K 

j 4>{x)z~'^Uz{x + ^)dx 
4>{x)uz{x)dx. 
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In this calculation we used the property of the Bloch solutions 

Uz{x + 7) = Z^fUzix). 

Therefore, 

< U,(p > = << fi{z),Uz >,(/)>, 

which concludes the proof of the theorem. □ 



4 Liouville-type theorem 

In this section we discuss Liouville theorems for periodic equations. We will 
consider at the moment an arbitrary linear elliptic operator P{x, D) with 
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smooth F-periodic coefficients which satisfies the assumptions made in Sec- 
tion ^ (as above, without loss of generahty we can reduce the consideration 
to the case T = Z"). 

Definition 22 We say that the Liouville theorem holds true for the operator 
P, if for any N gN the space Vi\i{P) of solutions of the equation Pu = in 
that can be estimated as 

||n|L2(^+^)<C(l + |7|)^for all 7 G T 

is finite dimensional. 

In the case when the Liouville theorem holds, we will be also interested 
in the dimensions of the spaces Vn^P) and in representations of their 
elements analogous to (|1.3|). 

The result below explains under what conditions on the operator P a 
Liouville-type theorem holds. These conditions will then be verified for some 
specific classes of operators. 

As was mentioned in the introduction, solutions representable as ( |1.3| ) are 
just Floquet solutions with zero quasimomentum. So, the Liouville theorem 
cited in the introduction states that any polynomially growing 
solution is a Floquet solution with a zero quasimomentum. Let us also 
mention that any Bloch solution e^^'^p{x) with a real quasimomentum k 
is automatically bounded. This means that the validity of the Liouville 
theorem for an operator P implies that the number of the real quasimomenta 
of solutions of the equation Pu = must be finite (modulo the action of 
the reciprocal lattice). In other words, the Fermi surface for P intersects 
the real space at a finite number of points (modulo the reciprocal lattice). 
In terms of the Floquet variety it means that the set := $p fl T is finite. 
We denote the cardinality of a set A by As the second statement of 

the next theorem shows, the finiteness of Z is in fact the only claim of the 
Liouville theorem. 

Theorem 23 1. The equation Pu = has a nonzero polynomially grow- 
ing solution if and only if it has a nonzero bounded Bloch solution, i.e. 
if and only if the intersection fl M" of the Fermi surface for P with 
the real space is not empty (or equivalently, Z = r\T ^ 'J)). 

2. The Liouville theorem holds for the operator P if and only if the inter- 
section Fp n M" is a finite set modulo the reciprocal lattice ( or equiv- 
alently, f^Z < 00). Moreover, if i^Z = 00 then the Liouville theorem 
does not hold even for bounded solutions, i.e., do = dimiVo) = 00. 
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3. 



If the Liouville theorem holds, then each solution u G Vn{P) can be 
represented as a finite sum of Floquet solutions: 



u{x) = E ^'pU^)- (4-1) 

(jeFpnR" lil<iv 

4. If the Liouville theorem holds, then for all N > we have 

dN < doqn,N < 00 , 

where qn,N is the dimension of the space of all polynomials of degree 
at most N in n variables. 

5. Assume that the Liouville theorem holds and that for each real quasi- 
momentum q (i.e., for each q G Fp OW^) the conditions of Theorem 
\Tdi are satisfied. Then for each N > the dimension d^ of the space 
Vn{P) is equal to the sum over q G {Fp n M")/r* of the dimensions 
of the spaces of Xq-harmonic polynomials (see Definition where Xq 
is the first nonzero homogeneous term in the Taylor expansion at the 
point q of the dispersion relation (band function) X{k). 

Proof: Statements 4 and 5 follow from 3 together with Lemma § and The- 
orem 10. So, we first prove statements 2 and 3 and conclude with the proof 
of the first statement. 

In order to prove 2 let us notice that if #Z = 00 then each point z = 
expifc G Z provides a bounded Bloch solution with the quasimomentum k, 
and these solutions are linearly independent. This means that the Liouville 
theorem cannot hold in this case. 

Assume now that < 00. We need to prove that the Liouville theorem 



and representation (4.1) hold true. Obviously, if u has a representation 
of the form (4J.), then u is of a polynomial growth. The proof that any 
polynomially growing solution is of the form ( |4.1| ) follows the same simple 
strategy as in the proofs of Theorem [l^ and as is in the proof of the main 
Floquet representatio n [30| , Theorem 3.2.1] (which, in turn, comes from the 
approach of |jl^ and |j9|]). As in the case with the "fundamental principle" 
(see [^] and |^^), it is more convenient to deal with a dual formulation, 
as it is done in |^0[. Namely, any polynomially growing solution u{x) can 
be interpreted in the dual way, as a functional on an appropriate functional 
space, which belongs to the cokernel of the dual operator P* . Consequently, 
a representation theorem for all such functionals must be obtained. In order 
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to make this idea precise, we need to introduce appropriate test functions 
spaces. 

Consider the Prechet spaces 

Cm = \^e Hl^Mll ^^vMln^iK+^y) (1 + ItI)"^ < oo, VAT 

Let the order of the operator P be m, then it is clear that P* maps con- 
tinuously Cjn into Cq. Due to the polynomial growth of u{x), the linear 
functional 

<u,(f)>= I u{x)(j){x)dx 

JRri 

is continuous on Cq. Since Pu = 0, one easily observes that u annihilates 
the range of the operator P* : Cm Co- We need now a Palcy- Wiener type 
theorem for the spaces Cm with respect to the Floquet transform. 

Lemma 24 1. The operator 

U : Cm ~^ C°°{T, £m) 

is an isomorphism, where C°°{T,£m) is the space of C°° sections of 
the bundle £m over the complex torus T, equipped with the standard 
topology. 

2. Under the transform U, the operator 

P* '■ Cm — > Co 

becomes the operator 1^(2) of multiplication by a holomorphic Fredholm 
morphism between the fiber bundles £m and Sq : 

C'^{T,£m)'^-^^ C^{T,£o). 

Here V{z) acts on each fiber of £m as the restriction to this fiber of 
the operator P* acting between H"^[K) and L'^{K). 

3. The operator V{z) is invertible for a point z if and only if ^ 

The next lemma is an analog of the classical theorem on the structure of 
distributions supported at a single point. Together with the previous lemma 
it essentially leads to the statement of the theorem. 
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Lemma 25 Let T be a C°° -manifold and V : T ^ L{Bi,B2) be a ex- 
junction with values in the space L{Bi,B2) of bounded linear operators be- 
tween Banach spaces Bi and i?2- Assume that for each z & T the operator 
V{z) is a Fredholm operator. Then 

1. IfV{z) is surjective for all points z in T, then the multiplication op- 
erator 



C°°{T, Bi 



Viz) 



C°^{T,B2) 



is surjective. 

2. If V{z) is surjective for all points z except a finite subset Z d T, 
then any continuous linear functional g on the space of smooth vector 
functions C^{T,B2) that annihilates the range of the multiplication 
operator 



C^{T, Bi) '^^ C°°(r,S2 



has the form 



<g,<P>=Y^ 

z€Z 



Dj,z{< 9j,z,(t>>) 



(4.2) 



Here gj,z are continuous linear functionals on B2, < gj,z,4> > denotes 
the duality between B2 and B2, Dj^z are linear differential operators 
on T, and iV G N. 

We are ready now to finish the proof of the nontrivial part of the third 



statement of Theorem 23 



If u is a solution of polynomial growth, it belongs, as it has been men- 
tioned already, to the cokernel of the operator P* : Cm Cq. After the 
Floquet transform we are dealing with the cokernel of the operator 

C^iT,£2)^-^^ C°°(T,fo). 



By Lemma 24, the only points z ^ T where V{z) is not invertible are those 
points where z~^ belongs to the Floquet variety. Since by our assumption 
the set Z = Tn<I> is finite, it follows that the operator function Viz) satisfies 



all the assumptions of Lemma 25. The fact that we are dealing with Banach 



bundles instead of fixed Banach spaces is irrelevant, since these bundles are 



28 



trivial. This means that we have the representation ( [4 .2] ) with gj G L^(T"). 
According to Lemma 0, functionals of the form (|4.2|) correspond under the 



inverse Floquet transform exactly to functions of the form (1.3). 

It remains to prove the first statement of the theorem. Let u be a 
polynomially growing solution. Assume that ^ = 0, i.e., the intersection of 
the Floquet variety ^ with the complex torus T is empty. Therefore, the last 



statement of Lemma 24 implies the invertibility of P{z) for all z E T. Now, 



the first statement of Lemma |2^ guarantees the surjectivity of the mapping 

and hence the absence of any nontrivial functionals on C°°{T,£q) that an- 
nihilate the image of this mapping. Since under the Floquet transform U, 
a polynomially growing solution u{x) is mapped to such a functional, we 
conclude that n = 0. □ 



Remark 26 The first statement of Theorem 23 is a part of the analog of 
the Bloch theorem provided in Theorem 4.3.1 of [|30| . Namely, the existence 
of a sub-exponentially (in particular, polynomially) growing solution implies 
the existence of a Bloch solution with a real quasimomentum, and hence the 
nonemptiness of the real Fermi variety. For completeness, we gave above an 
independent proof of this statement. 

One realizes now that the cases when a Liouville-type theorem holds in 
a nonvacuous way are extremely rare. Namely, Theorem ^ shows that this 
happens only when the Fermi variety touches the real subspace at a finite 
set of points (modulo the reciprocal lattice). This means in particular, that 
in the selfadjoint case, one should expect this to happen only at the edges 
of the spectral gaps. Although it is possible to imagine interior points of the 
spectrum where such a thing could occur, it is hard to believe that these 
cases could be anything more than accidents. 

One can expect the following conjecture to be true: 

Conjecture 27 Let P be a "generic" self-adjoint second order elliptic op- 
erator with periodic coefficients and (A_, A+) be a nontrivial gap in its spec- 
trum. Then each of the gap's endpoints is a unique (modulo the dual lattice) 
and nondegenerate extremum of a single band function \j[k). 

The validity of this conjecture together with Theorem ^ would imply 
that generically at the gap ends the dimension of the space V/v is equal to 
the dimension hn n of the space of all harmonic polynomials of order at 
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most in 71 variables. Unfortunately, the only known theorem of this kind 



is the recent result of |27], which states that generically a gap edge is an 
extremum of a single band function. 

At the bottom of the spectrum, however, much more is known. The 
theorem below combines some results of |17, 26, 41] with the statement of 
Theorem 23 to obtain the structure and dimension of the space of poly- 
nomially growing solutions in this case. Below the spectrum, the Liouville 
theorem holds vacuously, according to the first statement of Theorem 
and Theorem 5.5.1 in |30|. 



Theorem 28 1. Let H = —A + V{x) be a Schrddinger operator with a 
periodic real valued potential V E L^'/^(T"), r > n. Then the lowest 
hand function Xi{k) has a unique nondegenerate minimum Aq atk = 0. 
All other band functions are strictly greater than Aq. Every solution 



u G Vn{H — Aq) is representable in the form (1.5). The dimension of 
the space Vn{H — Aq) is equal to hn,N- 

2. Let V be like in the previous statement, then there exists e > such 
that for any periodic real valued magnetic potential A such that 

II^IIl'-(T") < ^ 

and 

A{x)dx = (4.3) 



the following statements hold true: The lowest hand function 
of the magnetic Schrddinger operator H = {iSI + A)"^ + V attains 
a unique nondegenerate minimum Aq at a point kQ. All other band 
functions are strictly greater than Aq. Every solution u € Vm{H — Aq) 
is representable in the Floquet form 

v{x) = e^'^o-^ xhAx) 

\j\<N 

with periodic functions Pj{x). The dimension of the space Vn (H — Aq) 
is equal to hn,N- 



3. Suppose that L is a second order elliptic operator of the form (2.5) 
such that Aq > 0. 
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/f A(0) =0 (i.e. G HJ, then the Liouville theorem holds and every 



solution u E Vm{L) is representable in the form The dimension 

of the space Vn{L) is equal to hn,N in the case when Aq = 0, and to 
Qn-i,N when Aq > 0. 

If A(0) > then the equation Lu = does not admit a nontrivial poly- 
nomially growing solution. So, the Liouville theorem holds vacuously. 



Proof: 1. The result of |26| says that the lowest band function Xi{k) has 
a unique nondegenerate minimum Aq at = and that all other band 
functions are strictly greater than Aq. Now Theorem implies the rest of 
the claims of this statement. 

2. When both the electric and magnetic potentials are sufficiently small, 
then the result of states that the lowest band function Ai(A:) of the 
magnetic Schrodinger operator H = (iV + A)'^ + V attains a unique non- 
degenerate minimum Aq at a point kQ, while all other band functions are 
strictly greater than Aq. This statement, however, can be easily extended 
to the case of arbitrary electric and small magnetic potential. Indeed, when 
the magnetic potential is equal to zero, one can refer, as in the previous 



case, to [26]. At this moment one has to use analyticity of the Bloch variety. 
Namely, the statement of Lemma Q (see also pO| , Theorem 4.4.2]) can be 
easily extended to include analyticity with respect to the potentials (see, for 
instance, [|l^). More precisely, there exists an entire function f{k,X,A,V) 
of all its arguments such that f{k, A, A,V) = is equivalent to 

{k, A) € i?(jv+A)2+y ) 

where Bh is the Bloch variety of the operator H. Now, the result of ]|2^ for 
A = together with the stated analyticity property imply the required fea- 
tures of the lowest band function for sufficiently small magnetic potentials. 



The last step is to use again Theorem 23. Note that the normalization 



always can be achieved by a gauge transformation which does not affect the 
spectrum and the Liouville property. 

3. The assumption A(0) > implies that the operator L admits a positive 
periodic supersolution. It follows from Lemma |l^ that the Fermi surface Fl 
can touch the real space only at the origin (modulo the reciprocal lattice 
r*) and in this case A(0) = 0. Therefore, by the first part of Theorem 
the Liouville Theorem holds vacuously if A(0) > 0. 

Suppose now that A(0) = 0. Lemma |l2| implies that if Aq > then the 
point = is a noncritical point of the dispersion relation, and if Aq = 
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then A; = is a nondegenerate extremum. Now Theorem |2^, as before, 
completes the proof. □ 

5 Proofs of the lemmas 

Proof of Lemma @: The first claim of the lemma corresponds to Theorem 



3.1.3 in |30|. In order to prove the second part of the lemma, let us fix a 
/cq £ C", and choose a closed subspace M C /^"^(T") complementary to the 
kernel of the operator P*{x,D — ko). Consider the (analytically depending 
on A; in a neighborhood of ko) subspace 

n(A:) := P*{x, D - k){M) C ^^(T"). 

and 

M := mko)]^ . 

Then dim{J\f) = a^Q, and for values of k close to /cq the space Af remains a 
complementary subspace to Il{k). Representing the operators P*{x, D — k) 
in the matrix form according to the decompositions 

fjm^jn^ = M®Ker P*{x, D - ko) 

and 



L2(T") = n(A:)eAr, 



we get the matrix 



B{k) * 
C{k) 

where B{k) is an invertible analytic operator function, and C{k) is an ana- 
lytic matrix function of the size o^q x a^^. Here a^^ is the dimension of the 
kernel of the operator P*{x,D — ko). (Notice that ak^ = if indP = 0, 
which is true for instance, when dealing with scalar elliptic operators, due 
to the Atiyah-Singer theorem.) Now, the space of all distributions orthog- 
onal to the range of P* and supported at exp(— i/co) reduces to the space 
of all distributions supported at ko, acting on C^* -valued vector functions, 
and orthogonal to the range of the operator of multiplication by C{k). If 
we drop the orthogonality condition, the dimension of the space of all such 
distributions of order at most N is obviously equal to akqn,N, which proves 
the estimate. We point out that a direct proof of this estimate for scalar 
operators can be also easily derived using the Leibnitz's rule. □ 
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Proof of Lemma 12: Statements 1 through 3 of the lemma are contained 
in 1^^, except the statement that the geometric rather than the algebraic 
multiplicity of the eigenvalue A(^) is equal to one. The latter follows easily 



from Lemma 5.2 of |3C]. Alternatively, it can be deduced from general 



theorems on positive operators defined on an ordered Banach space (see for 



instance, |29, Theorem 2.10]). Statement 4 is proven in [41, Theorem 5]. □ 



Proof of Lemma 13: Statements 1-3 follow from the results of |2|, 



while statements 4-5 follow from Theorem 5]. □ 
Proof of Lemma |14| : Consider the following family of operators on the 



torus: L{x,D — i^) — A(^). It follows from Lemma 12 that this family 
is analytic in a complex neighborhood W of the set G and its values are 
Fredholm operators between the appropriate Sobolev spaces. The same 
lemma implies that the dimension of the kernel of all these operators is 
equal to 1. Hence, these kernels form an analytic fiber bundle over W (see 
Theorem 1.6.13 and the corresponding references in |^^). One can always 
assume that the domain W is convex (in the geometric sense). Then the 
kernel bundle (as all vector bundles on W) is topologically trivial. Since 
W, being convex, is a domain of holomorphy (see for instance Corollary 



2.5.6 in |24|), therefore, the result of (an instance of the so called Oka's 
principle) implies that the bundle is also analytically trivial. This means the 
existence of a nowhere zero analytic section u^. Positivity of for G H 
can be achieved as follows. Let us choose any nonzero analytic solution as 
above. Then for some small neighborhood Wi cWofG, we have n^(0) ^ 0. 
So, we may normalize by dividing it by u^{0). The resulting solution is 
clearly positive for ^ € H. □ 

Proof of Lemma |15| : 1. Let u{x) = e^^'^p{x) be a nonzero Bloch solution, 
where p{x) is a F-periodic function, and k E F f] T. Assume first that 

o 

Im k E — G, so, Aq > 0. We need to prove that Re u = Im u = 0. We show 
for instance, that ui := Reu = 0. Suppose that ui ^ 0. We may assume 
that ni(xi) > 0, for some xi € M". Consider the positive solution 

v{x) = / Ui:{x)da{S,), 



where da is the (n — l)-dimensional surface area element on H. For every 
M > there exists R > such that v{x) - Mui{x) > for all \x\ > R. 
By the generalized maximum principle, v{x) > Mui{x) in M". Since M is 
arbitrarily large and ui{xi) > 0, we arrived at a contradiction. Note that 
this argument applies also to any Floquet solution with a quasimomentum 
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k such that Im k — G- 

Suppose now that Im k G — H and A > 0. Clearly, it is enough to show 
that there exists a real constant C and ^ G H such that ui := Reu = Cu^. 
Let ^ = —Imk. Then for a sufficiently small e > the function Ve '■= 
^ — eui is a positive solution of the equation Lu = 0, which is smaller 
than u^. Recall that is a minimal positive solution of the equation 
Lu = 0. Therefore, there exists c > such that = cu^, which implies 
that ni = Cu^ for some (7 € M. 

2. Consider the zero set Fi of the analytic function A(ik) in a small complex 
neighborhood of — iH. Since Aq > 0, it follows that the gradient of A{ik) is 
not zero on —i'E. Therefore, Fi is a smooth analytic variety. We will show 
that the Fermi surface F coincides with Fi in a neighborhood of —i'E, which 
will conclude the proof of the lemma. Indeed, obviously Fi C F. Consider 
a point ko = —i^o G —iE. By Lemma |l2|, zero is a simple eigenvalue of the 
operator L(x, D+ko) = L{x, D—i^o). This means that the spectral projector 
that corresponds to a neighborhood of zero is one-dimensional for all complex 
k close to /cq- We conclude that for all A; in a complex neighborhood of — 
there is exactly one eigenvalue close to zero of the operator L{x,D + k). By 
Lemma ^3], we know this eigenvalue, namely K{ik). Let now k belongs to a 
small neighborhood of —i'E and assume that k ^ Fi. Then K{ik) ^ 0, and 
hence zero cannot be the eigenvalue of D + h). This means that k does 
not belong to the Fermi surface F. □ 

Proof of Lemma |17| : 1. If c ^ 0, the assertion of the lemma follows from 
Theorem 4.5]. On the other hand, if c = 0, then G — S, and in 



particular, G {—G). It follows from Lemma 15 that any Bloch solution 



with a real quasimomentum is the constant solution. 

2. This assertion follows directly from the part 1 using the operator ijj^'^Lilj. 

□ 

Proof of Lemm a |20| : The second statement of the lemma coincides with 
Theorem 2.2.3 in [Q. So, we need to prove only the first statement. 

Let if G Wm,e- We will show that the series ( |2.3[ ) converges uniformly on 
compacta in V* as a series of functions on V* with values in H'^{K). This 
would imply that Uip G TiV*, H"^{K)), and that the corresponding (one-to- 
one) mapping W : Wm,e T{V* , H"^{K)) is continuous. Let < 6 < 6i < e. 
Let z = expik G which means that, Imk & Gs- We have 



34 



I •) 1 1^.(^) < ^ 1 1^1 e-^-^-^ = 1 1^1 e^- '^•^ 

7Gr 

We need to check now that the mapping U acts from Wm,£ into T{y*,£rn)- 
This amounts to showing that satisfies the appropriate Floquet boundary 
conditions and hence is in fact a section of the sub-bundle £m C V^^ x (K) . 
This is a straightforward calculation (see also Theorem 2.2.2 in |30|] ). 

On the other hand, let us assume that s{z) € T{V* ,£m)- If 2; = expik, 
then s as a function of k is periodic with respect to the reciprocal lattice T*. 
Expanding it into the Fourier series, we get 

7Gr 

where s-y € H'^{K). We can now define a function (p on M" such that 
(p{x — 7) = 5^(2;) for X € and 7 G F. 

The function (p belongs to H"^ in the interior of each of the cubes K + 
One only needs to check that it belongs to H"^^ at the boundary points of 
these cubes. The requirement that s{z) is a section of the bundle Em rather 
than just of the bundle V* x H^{K) does exactly this (see the discussion 
at the top of page 96 in [pO|]). 

It remains to show that ip £ Wm.e- We use the standard formulas for the 
Fourier coefficients to get 

^{.-j)=s, = j^ I .(e^('5+-))e-W«)-7rf^, VaGG,, 

B 

where B is the first Brillouin zone, and we write z = exp ik = exp{i{P+ia)). 
Note that 

IMn'^iK+y) ^ ™a? \ \s{z)\\Hm^K) e"'^"^'' Va G , (5.1) 
and therefore, 

M\H-iK+y) < max \\siz)\\^^^j,^e-(^(^/\^\^+'^K (5.2) 
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This implies immediately that 



< ^ >„,5= sup { 1 |^™(^+^) e('^(^/l^l)+^)l^l } (5.3) 



7Gr 

< C max ||s(z)||H-™,r^^ supe"*-'^^^'^^''''' < oo, 



if 5i >(^. □ 



Proof of Lemma 21: The statement of this lemma is established in a 
much more general situation at the beginning of the proof of Theorem 1.7.1 
in However, for the sake of completeness we provide here the proof 

for our simpler particular situation. First of all, the sequence of the lemma 
is a complex (i.e., the composition of any two consecutive operators in it 
is equal to zero). One needs to prove this only in the second term of the 



sequence, where it follows immediately from the equality (3^). Indeed, since 



Uz solves the equation Lu = 0, (3^) followed by integration by parts proves 
the statement. 

Let us turn to the exactness. We need to prove it in the second and third 
terms of the sequence. Consider the second term. Let f{z, x) G TiV* ,£q) be 
such that tf{z) = 0. This means that for any z G the function f{z~^, •) 
is orthogonal to the Bloch solution of the equation Lu = 0. We need to 
show that g{z) = C{z~^)~^ f (z) is analytic, which will mean that / belongs 
to the range of C. The function g{z~^) is automatically analytic outside of 
<I>e, so we only need to make sure that it does not develop any singularities 
at this subset. We will show that all the necessary and sufficient conditions 
for the analyticity of g have the form of orthogonality of values of / at 
certain points to certain functionals. This would resolve the issue, since all 
such possible orthogonality conditions are the orthogonality of f{z~^) to the 
kernel of L on Bloch functions with a quasimomentum z, and hence to the 
vanishing of tf{z). As it was shown in the proof of Theorem 3.3.1 in 
pages 113-114], the inverse operator to C{z~^) is the ratio of two analytic 
functions: 

C{z-')-' = B{z)/A{z), 

where B{z) is an analytic function with values in bounded operators from 
L2(T"') to ff^(T"), and A(z) is a scalar analytic function, which is a regular- 
ized determinant of £{z~^)C{zq^)~^ for some point zq where the operator 
is invertible. Such regularized determinants are determined in the stan- 
dard way by the eigenvalues of the corresponding operators (see for instance 



Section 2 of Chapter IV in [20| for general definitions and properties of regu- 



larized determinants, and for our particular situation the proof of Theorem 
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3.1.7 and related discussion in Section 1.2 in [^]). The simplicity of the 
eigenvalue A(^) (Lemma ^) implies that if we introduce instead of z the 
coordinate ^ such that z = exp^, then A(exp^) = A(^)Ai(^), where Ai(^) 
is an analytic function with no zeros in the domain under our consideration. 
We recall now that A has simple zeros. Hence, the necessary and sufficient 
condition for / to belong to the range of the operator C on the space of ana- 
lytic sections is that the vector- function B{z)f{z) vanishes on the set of the 
zeros of A. These conditions obviously have the form of the orthogonality 
of values of / to some functionals. As it was explained above, this implies 
exactness at the second term of the sequence. 

Let us turn now to proving the exactness at the third term. We need 
to show that arbitrary analytic function on can be obtained as tf{z) for 
some / G T{V;,£,). 

Let us denote by the manifold 

q>* = {z\z-^ € «>J. 
Consider the restriction mapping 

r(v;*,^o)^r($:,^o). (5.4) 

Notice that <I>* is an analytic subset in V* and that and V* are domains 
of holomorphy. The latter can be easily proven using power test functions z"" 
with integer (but not necessarily nonnegative) powers a (a similar derivation 
can be found in the proof of the implication {iii) — > (i) of Corollary 2.5.8 
in 1^^). Then Corollary 1 of the Bishop's theorem p3| . Theorem 3.3] (see 
the original theorem in |[l^) claims that the restriction mapping ( |5.4D is 
surjective (recall that the bundle Eq is trivial). Hence, it is sufficient to 
prove that the mapping 

t:r($„^o)^w(^.). 

defined as 

t/(z) =< /(z, ■),Uz >= j f{z,x)uz{x)dx 

is surjective. Consider the continuous operator T{z) : Lp'{K) C defined 
as T{z)y =< y,Uz >= / y{x)uz{x)dx. Since Uz is not zero, this operator is 

surjective. It is clear that it depends analytically on z. According to Allan's 
theorem (see |^] or Theorem 4.4 in |4^]), since $e is a Stein manifold, there 
exists an analytic right inverse operator R[z). Now, given ^(z) £ 7^(<I>e), 
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the function g{z) = R{z)(p{z) satisfies tg = (p. Tliis proves the surjectivity 
that we need. 

The last statement of the lemma about the mapping t being a topological 
homomorphism is just the open mapping theorem. □ 



Proof of Lemma 24: 1. We first show that the operator U maps con- 
tinuously the space Cm into C°° {T , H™- {K)) . Indeed, if E Cm-, then 
ll¥'lli?™(_K'+7) decays faster than any power of |7|. This together with (|2.3| ) 
leads to the immediate conclusion that Uip belongs to C°°{T,H"^{K)) and 
to the continuity of the corresponding mapping. Since U(p is a section of 
the sub-bundle (see the Section 2.2 in [Q), this gives us the needed 
conclusion. Conversely, let 

s{z) E C7°°(r,^™) C C°°{T,H"'{K)). 

One can expand the -ff™'(i^')-valued function s{z) into the Fourier series: 

s{z) = s.yz"', z 

Here s-y G H^{K). Standard estimates of the Fourier coefficients of smooth 
functions apply, which show that decays faster than any power of |7|. 
Let us define now a function (j) on such that — 7) = s^{x) for x £ K 
and 7 G r. The additional information that s is a section of the sub-bundle 
leads (as in Q, page 96]) to the conclusion that G Hl^^iR""). This 
implies that <j) G Cm and finishes the proof of the first statement of the 
lemma. 

Statements 2 and 3 are correspondingly parts of Theorem 2.2.3 and 3.1.5 



ofl3Cl. □ 



Proof of Lemma |25f The first statement is rather obvious. Indeed, the 
statement is local, and locally one can construct a smooth one-sided in- 
verse. The second statement can be proven like the similar statement in 
Corollary 1.7.2]. For completeness, we provide the scheme of the proof here. 
Under the conditions of the second statement of the lemma, it is easy to see 
that any functional annihilating the range of the operator of multiplication 
by V{z) must be supported at the finite set Z where V{z) is not surjec- 
tive. This also reduces the considerations to a neighborhood U of a point 
zq G Z. Using the Fredholm property, one can find a closed subspace M of 
finite codimension in Bi such that the operators Viz) have zero kernel on 
M for all z £ U (see the corresponding lemma in Q, or Lemma 1.2.11 and 
Remark 2 below it in [^]). Now the problem reduces to a similar one on 
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a finite-dimensional space, where a standard representation of distributions 
supported at a point implies (|4.2|) . □ 



6 Further remarks 

Remarks 6.1 1. Throughout the paper, we have assumed for simplicity 
that all the coefficients of the operators P and P* are C°°-smooth. In 
fact, we do not need such a restrictive assumption (see the discussion in 
Section 3.4.D]). For example, a sufficient (but not necessary) condition for all 
the statements of Section |3| to hold true is that the coefficients of L and L* are 
Holder continuous. Actually, even less is needed. For instance, conditions 
imposed on the Schrodinger operators in Theorem [2^ are sufficient. It is clear 
that the conditions on the coefficients could be significantly relaxed, if the 
operators were considered in the weak sense, or by means of their quadratic 
forms. This should not change the general techniques of the proofs. We did 
not intend, however, to find the optimal requirements on the coefficients for 
all our results to hold. 

2. It should be possible to describe the class of solutions of the equation Lu = 
that are representable by a distribution rather than by a hyper function. 
We plan to address this problem elsewhere. 

3. The Liouville theorem can probably be extended to systems of equations 
(for instance, to the Maxwell system). In this case one would face the 
problems of a possibly nonzero index of the corresponding operator and of 
multiple eigenvalues (the latter can also occur for scalar operators). We 
believe that the technique of this paper might be adjusted to handle some 



of these situations. The extensions of the result of |26] to the Pauli and 
Maxwell operators obtained in |0] and |0] would provide examples where 
the needed information on the behavior of the dispersion relations at the 
bottom of the spectrum is available. 



A Appendix 



In this appendix we present an alternative proof of the third statement of 
Theorem 28 in the case when either Aq 

< iV < 

theory [^5| and imitates the proof of Theorem 2 in |38], where L is assumed 



and iV > 0, or Aq > and 
The proof relies on some basic notions of homogenization 



to be an operator in divergence form. Therefore, we skip some details which 



are essentially the same as in |38 
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We need to recall some basic definitions from homogenization theory (see, 
for example, |10, ^). Suppose that L is a second order elliptic operator of 
the form 



L = - ^ aij{x)didj + ^ bi{x)di , (A.l) 

i,j=l 1=1 

with periodic coefficients and denote the positive matrix {aij{x)} by A{x) 
and the periodic vector (61, . . . , &n)^ by b. Let tp be the positive normalized 
periodic solution of the equation L*u = 0. Let ^(x) = (^i(x) . . . , ^'^(x))^ 
be a solution of the equation 

L^ = -b{x)+[ b{x)^{x)dx inT". (A.2) 

Consider the matrix 

Q = ■■= f (1 + V^fA{x){I + V^)^{x) dx , (A.3) 

were / is the identity matrix. The operator Q := —Yl'ij=i1ij^i^j called 
the homogenized operator of the operator L, and the positive matrix Q = 
{qij} is called the homogenized matrix (see, p5| . Section 2.5]). 

The following lemma, which is actually a new formulation of |41, Theo- 



rem 5]), establishes a connection between the function A and homogenization 
theory. 



Lemma A.l Let L be an operator of the form ^2.f\ ) and suppose that ^ G H. 
Let and itl^ be the positive Bloch solutions of the equations Lu = and 
L*u = 0, respectively. Denote by tp the periodic function u^u*_^. Consider 
the operator 

n n 

L = {u^(x)y^ Lu^{x) = - ^ aij{x)didj + ^ bi{x)di , (A.4) 

i,j=l i=l 

let 

n 

Q = - ^ Qijdidj (A.5) 

be the homogenized operator of the operator L, and Q = {qij} be the homog- 
enized matrix. Then ip is the principal eigenfunction of the operator L* on 
the torus T" with an eigenvalue 0. Moreover, Hess(A(^)) = — Q. 
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Proof: The first statement of the lemma can be checked easily while the 
second statement follows directly from the formula in [pT| , Theorem 5], and 
the definition of the homogenized operator. □ 



Proof of a part of the third statement of Theorem g8 

may assume that LI = 0, so, 



We clearly 



L = - ^ aij{x)didj + y^^bi{x)di. 



i,j=l i=l 

We denote by ^ the normalized positive solution of the equation L*u = 
in ¥"■. Let ^ be a solution of the system (A. 2), and Q be the homogenized 
operator of the operator L. 

Assume first that Aq > 0. The case = is trivial, and follows from 
Theorem 23 and Lemma 0. Let = 1. Recall that according to Theorem 



23, di < n + 1. Moreover, by Theorem 23 and the Leibnitz's rule, a (real) 
solution of linear growth is of the form 



u{x) = ''^^CLjXj + (p{x), 



where aj € 



and (f) is periodic. 
By lemma 13, Aq = if and only if for every 1 < k < n 



a 



hj{x)'4}{x) dx = . 



(A.6) 



For 1 < J < n, we write an "Ansatz" for a solution of linear growth of the 
form 



Fj{x) = Xj + (pjix), 



(A.7) 



where (f)j is a periodic function. Clearly, Fj is a solution of Lu = in 
if and only if (pjix) solves the nonhomogeneous equation Lu 
By the Fredholm alternative, this equation is solvable in T" if and only if 
'.j = which holds true for all 1 < j < n, if and only if Aq = (and in this 



bj in T". 



a 



case, 



^'j, see ([A^). Therefore, di = n + 1 if Aq = 0, and di < n + 1 



if Ao > 0. 

In order to finish the proof for = 1, we need to prove that if Aq > 0, 
then di > n. Without loss of generality, we may assume that a„ 7^ 0. We 
construct (n — 1) linearly independent solutions of linear growth of the form 



Fj{x) = Xj - aj{an) x„ + (j)j{x), 
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where 1 < j < n — 1, and (pj solves the equation Lu = —bj+ aj{an)~^bn- 
Note that these (n — 1) equations are solvable and therefore, di > n. 

For N > 2, we assume that Aq = 0. Recall that if u E V/v then by 
Theorem 23 and the Leibnitz's rule 

u{x) = u^^\x) + x^pi^{x), 

\u\<N 

where 

\u\=N 

and pi, are periodic functions if < N, and pi, E M, if = A^. 
Claim: Assume that Aq = 0. Then for all N >0 

Q'u(^)=0. (A.8) 

In particular, cIn < hn,N- 

Proof of the claim: Assume first that N = 2. Then u G V2 is of the form 

1 " 
u{x) = -{Cx ■ x) + YxjPj{x) +po{x), 

i=i 

where C is a constant symmetric matrix, and pQ,pi,... ,pn are periodic 
functions. 

A direct calculation shows that the vector p = {pi, . . . ,pn)'^ must satisfy 
the equation Lp = —Cb which is solvable since Aq = 0. Therefore, p = 
(up to a constant vector). Also, po must satisfy 

Lpo = f := tr{A{I + 2V^^)C7^) - b ■ . 

The compatibility condition for this equation is / f{x)tlj{x) dx = which 

T" 

after some calculations implies that 

tr(QC7^) = 0, 

where Q is the homogenized matrix of the operator L (see ( |A.3| )). Since 
u^"^^ := ^{Cx ■ x) is a homogeneous polynomial of degree 2, it follows that 
Qii^^) = tr{QC'^). Therefore, u^"^^ solves the equation Qu = 0. Thus, the 
case = 2 is settled. 

For > 2, we proceed by induction as in |^^. Namely, assume that 
the claim (|A.8D has been proven for — 1, and let u S Vat. Let Aj be the 
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difference operator Aif{x) := f{x + Cj) — f{x), where is tfie i-th. vector 

of the standard basis of M", and 1 < i < n. Then Vi := AjU € Vat-i and 
the leading part of Vi is given by {Aiu)^^~^^ = diU^^\ By the induction 
hypothesis, (5((Ajti)(^^^)) = 0. Therefore, 

5,(Qu(^)) = Q(ainW) = Q((A,ti)(^-i)) = l<i<n. 

Hence, Qu^^^ = const. , and since Qu^^^ is homogeneous of degree iV — 2 > 
0, we obtain that Qu^^^ = 0, and the claim is proved. 

It remains to prove that d]\f > hn^N- So, for any homogeneous polynomial 
h of degree N which is Q-harmonic, we need to find a solution u G Vn such 
that u^^^ = h. Let u G Vn and e > 0. Consider the function 

\u\<N 

which tends to w/^-* as £ ^ 0. We consider x and y = | as independent 
variables and write 

U{x,y,e):= ^ s^-^^^x^p^iy) = Uo{x) + sUi{x,y) + ■ ■ ■ + e''UN{x,y) . 

Then the equation L{x,dx)u = implies that 

{Lo + eLi+e'^L2)U = 0, 

where 

n n n 

Lo = L{y,dy) ; Li = -2 ^ aij(y)a^. +^ 6i(y)5x, ; L2 = -Yaij{y)dl.^^.. 

i,j=l 1=1 i,j=l 

We look for a formal differential operator 

oo 

3=0 V 

where <?!>;/ (y) are periodic functions and 0o = 1- This operator should satisfy 

n 

(Lo + eLi + e^L^)^ = M + Lo(y, dy) - 2£ ^ aij{y)dx,dy. , (A.9) 

2j = l 
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where the formal operator 

oo 

3=2 \u\>2 

has constant coefficients. 

Comparing the coefficients of in ( [A .91 ) yields the following equations 
(the equation for s = is automatically satisfied). 

n 

Lo$i +Li = -2Y^ aij{y)d^^dy^ , s = 1 , (A. 10) 

Lo^s + Li^s-i + L2^s-2 = Ms, s>2. (A.ll) 



It is easily checked that for s = 1 the functions (pjiy) of Equation 
are the corresponding solutions for <I>i. Also, Equation ( |A.11| ) for s = 2 is 
solvable if M2 = Q, where Q is the homogenized operator of L. Similarly, 
the constant coefficients of the operator Mg , s > 2, are determined by the 
compatibility condition for Equation (A.ll) with s > 2. 

Let R : "P — > P be a linear right inverse of the homogenized opera- 
tor Q that preserves the homogeneity of polynomials. Consider the formal 
operator A which is defined by the equation 

00 

A-I = RY,£^Mj+2, 

and let be its unique formal inverse. Note that e^M2A = M. 

Let Uq{x) be a given homogeneous polynomial of degree which solves 
the equation Qu = 0, and let V{x) := A~^Uo{x). We have 

MV = e^M2AV = e^AhUo = 0. 

Define U{x,y,e) := ^A^^Uq = ^V, and denote u{x) := U{x,x,l). By 
inspection, u has a polynomial growth of order A^, and u^^\x) = Uo{x). 
Moreover, 

(Lo + eLi + e^L2)U = (Lq + eLi + e^L2)<^V 

n 

= MV + Lo{y, dy)V{x) + 2e aij{y)d,^dyy{x) = , 
and ^A~^ is the desired mapping. □ 
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Remark A. 2 1. Let Fj be the solutions of hnear growth defined by Equa- 



tion ( A. 7 ). A. Ancona g proved that the map F{x) = . . . ,F„(x)) is 

a diffeomorphism on if n < 2, while for n > 2 this map is not necessarily 
a diffeomorphism. 

2. Assume that LI = and Aq = 0. Let A(^) = ^|,^|>2 0^C'^ be the Taylor 
expansion of the function A. We conjecture that = m^, where rrii, are 
the coefficients of the operator M. 
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